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The statistics of 2-dimensional turbulence exhibit a riddle: the scaling exponents in the regime 
of inverse energy cascade agree with the K41 theory of turbulence far from equilibrium, but the 
probability distribution functions are close to Gaussian like in equilibrium. The skewness <S = 
Ss(R) I S^ 2 (R) was measured as <S cxp ~ 0.03. This contradiction is lifted by understanding that 
2-dimensional turbulence is not far from a situation with equi-partition of enstrophy, which exist 
as true thermodynamic equilibrium with K41 exponents in space dimension of d — | . We evaluate 
theoretically the skewness S(d) in dimensions | < d < 2, show that S(d) = at d = |, and that it 
remains as small as <S exp in 2-dimensions. 
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Two dimensional (2D) turbulence is not realized in na- 
ture or the laboratory, but only in computer simulations 
. Nevertheless it serves as an idealized model for a va- 
riety of natural flow phenomena, like geophysical flows in 
the atmosphere, oceans and magnetosphere. Experimen- 
tal setups that are close to 2D turbulence were realized 
in a number of laboratories § , and the advent of faster 
computers allows precise simulations of the 2D Navier- 
Stokes equations ||. In this Letter we are interested 
in the statistical characteristics of 2D turbulence, in the 
probability distribution functions and spectra of velocity 
differences. The velocity difference across a scale R is 
written in terms of the velocity field u(r,t), 

w(r,R,t) = u(r + R,t)-u{r,t) . (1) 

Usually one measures the longitudinal component, 

w e {r, R, t) = [u(r + R,t) - u{r, t)] ■ R/R , (2) 

the probability distribution function (pdf) of this object, 
denoted as P[we(r, R)], as well as moments of this pdf, 
like the second and third order structure functions 

S 2 (R) = (w 2 {r,R,t)) , S 3 (R) = (w 3 e {r,R,t)) . (3) 

Here the average is over space and time. In stationary 
homogeneous and isotropic ensembles the structure func- 
tion depends on R only. 

It is well known that in 2D turbulence the properties of 
these objects depends on whether R is larger or smaller 
than the scale L at which energy is injected into the sys- 
tem. For R^$> L (but smaller than the outer boundaries 
of the system) one observes an inverse cascade of energy, 
and S2(R) is found to scale with an exponent in agree- 
ment with the Kolmogorov 1941 theory, i.e. 

S 2 (R) = C 2 (eR) 2 ' 3 , i?>L. (4) 

Here C 2 is a dimensionless coefficient of the order of unity, 
and e is the mean energy flux per unit time and mass. 
For R <C L (but larger than the dissipative scale) one 
observes a direct cascade of enstrophy, with an exponent 



close to 2, but with some logarithmic corrections. Recent 
experiments and simulations lend strong support to Eq. 
([|); on the face of it this indicates that the system is 
very far from equilibrium, where cquipartition of energy 
is expected. On the other hand, the same experiments 
and simulations indicate that P[wi(r, R)] appears almost 
Gaussian, as if the system were very close to equilibrium. 
Quantitatively one measures the skewness 

S = S 3 (R)/S 3 2 /2 (R) , (5) 

with the results that <5> cxp « 0.03 <C 1. This seeming 
contradiction and its resolution are the subjects of this 
Letter. 

The basic idea of this Letter is that there exists one 
space dimension for which these observations are not at 
all in contradiction. This is d — |, in which there ex- 
ists an equilibrium state with equipartition of enstrophy, 
where the scaling expected for 52 (R) is exactly (Q), and 
where the skewness S — 0. By examining turbulence in 
|+e dimensions, we establish that throughout the range 
< e < | the situation remains very close to the one 
seen at e = ; the value of S(d) provides a natural small 
parameter to characterize the distance from d = | . This 
parameter remains very small up to 2-dimensions. We 
thus interpret the statistics in the inverse cascade regime 
of 2D turbulence as a state very close to equilibrium. 

It should be said at this point that attempts to con- 
nect turbulence to equilibrium statistical mechanics were 
made before. Well known are the theories advanced by 
Onsager Q, Hopf [|| and Lee ||, and see ref. for a 
review. There ideas from equilibrium statistical mechan- 
ics were proposed to explain turbulence; as summarized 
in 0| , the consensus is that such ideas were not relevant 
for forced turbulence, but may be relevant for the final 
stages of the temporal evolution of decaying turbulence. 

Next there were attempts to connect the physics in 2- 
dimensions to special properties of other, non physical 
dimensions. Of particular influence were ideas advanced 
by Fournier and Frisch 0, who identified d w 2.05 as 
the dimension at which the direction of the energy flux 
changes sign. It was hoped that this may provide a con- 
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venient starting point for perturbative expansions in the 
magnitude of the flux. A similar point was identified in 
shell models of turbulence, as a function of a parame- 
ter ||. It turned out however that the fluctuations were 
not small at that point, the flux changed sign discon- 
tinuously, and perturbative theories did not yield useful 
insights. In our own work on shell models M we iden- 
tified another point in parameter space where a quasi- 
equilibrium state with equi-partition of "enstrophy" co- 
incided with an energy spectrum in agreement with Kol- 
mogorov scaling exponents as in (g). This has led us 
to seek an 1-parameter lifting of the Navier-Stokes equa- 
tions where a similar phenomenon exists. We were thus 
led to examining turbulence in d = | and its vicinity. 
We argue below that there is a fundamental difference 
between the properties of turbulence in d — | (where the 
flux vanishes) and d ss 2.05 (where the flux discontinu- 
ously changes sign). At the former point the statistics 
is exactly Gaussian, with equipartition of the enstrophy. 
The latter has no such distinction, and the statistics may 
differ strongly from Gaussian. 

We begin by generalizing 2D turbulence to d < 2 di- 
mensions. In doing so we preserve the d-dimensional en- 
ergy and enstrophy, in distinction with other such exten- 
sions that fail to do so. The Navier-Stokes equations for 
the velocity field u(r,t) read 

du/dt + (u ■ V)« - vV 2 u + Vp= f, V • u = , (6) 

where v is the kinematic viscosity, p is the pressure, and f 
is a force which maintains the flow. In d = 2 it is natural 
to consider the vorticity field u>(r,t) = V x u(r,t). For 
v = / = the curl of Eq. (||) is du:/dt = V x (u x u). 
Introduce the right-handed coordinate system (x\ = x, 
x 2 = y and £3 = z) in which for 2D flows U3 = and 
u(r,t) = u(x,t), where x — {xx,x<z}. The vorticity has 
the single component W3 for which Eq. (pj) reduces to 



du}/dt+(u-V)uj = Q 



(7) 



This equation has two quadratic integral of motions, the 
(kinetic) energy E and the enstrophy H: 

E = i/" u 2 (x 7 t)d 2 x, H =\j u 2 {x 7 t)d 2 x . (8) 

The velocity and vorticity of a 2D flow may be derived 
from the stream function tp(x,t): 

u(x,t) = -V X ztp{x,t) , u(x, t) = -V 2 ijj(x, t) , (9) 

where z is a unit vector orthogonal to the x plane, 
and V 2 is the Laplacian operator in the plane. In fc- 
representation, a(k, t) = k J dr exp[—i(r-k))ip(r,t) . The 
Fourier transforms of u(x, t) and of u(x, t), respectively 
u(k,t) and u)(k,t), are expressed in terms a(k,t): 

u(k,t) = i(z x fe) a(k,t) , w(k,t) = -ka(k,t) , (10) 



where k = k/k. Now, by Eqs. (0) and @) 
da(k,t) _ 1 fd 2 qd 2 p 



dt 2 J (2tt) 2 6 ( k + 1 + P) V kqp a *(q,t) a *(P,t), 
Vkqp = S kqp {p 2 - q 2 )/2kqp, S kqp = 2qp smcf> vq , (11) 

Skqp — $pkq — Sqpfc — $kpq — Sqkp — Spqfc , 

|S*p,| - ^(fcV + q 2 P 2 + P 2 k 2 ) -tf-qi-p* . 

Here the interaction amplitude (or "vertex") Vk qp is ex- 
pressed via Skqp] \Skqp\/4: is the area of the triangle 
formed by the vectors fc , q and p. <j} pq = <j) p — (f> q ; (j>k , 4> q 
and cj) p arc the angles in the triangle plane between the 
xi-axis and the vectors fc, q and p respectively. 
The vertex Vk qp satisfies two Jacoby identities Jl|] 



{Vkqp + Vpkq + Vqpk) 
{k 2 V k qp + p 2 Vp k q + q 2 V qp k) 



(12) 
(13) 



The first one guarantees the conservation of energy in 
the inviscid forceless limit, while the second conserves 
the enstrophy. In terms of a(k,t) Eqs. M) read 



\a(k,t)\ 2 d 2 k 



k 2 \a{k,t)\ 2 d 2 k 
2(>p 



(14) 



For a statistical description in d-dimensions we con- 
sider the 2nd and 3rd order correlation functions of a. a* : 



(2ir) d 5(k + q)n k (t) 
(2TT) d 5(k + q+p)F k qp(t) 



(a(k,t)a(q,t)) , (15) 
(a(k,t)a(q,t)a(p,t)) . (16) 



In isotropic systems, we do not need to carry the bold- 
face fc index in n k and in Ff- qp . In terms of rik we define 
the volume densities in d dimensions: 



d d k - N= H_ 
2(2ir) d nkl n ~ V ~ 



d d k 
2(2tt)' 



k 2 n k ■ 



(17) 

Since rik and k 2 nk serve as the energy and enstrophy den- 
sities respectively in fc-space, the thermodynamic equilib- 
rium can be achieved with equi-partition in any of these 
quantities: 



,e0 



JiO 



A, 



energy equi-partition; 



(18) 



Ah/k 2 , enstrophy equi-partition 



In such a state of thermodynamic equilibrium unavoid- 
ably all fluxes vanish and the pdf of the velocity differ- 
ences is Gaussian. Usually in the theory of turbulence 
one rather considers flux equilibria, which in the present 
situation can be with energy flux or enstrophy flux. Di- 
mensional analysis for such flux equilibria predicts 

n% = C e (d)e 2/3 k~ x ? , a: e = d+f, energy flux; (19) 
n k l = C h {d)h 2/z k~ Xh , x h = d + 2, enstrophy flux(20) 

Here h is the mean enstrophy flux per unit time and mass, 
whereas C s {d) and Ch{d) are d-dependent dimensionless 
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coefficients. In terms of Sa(-R) these results are in agree- 
ment with (Q) for all d for energy flux equilibrium. For 
enstrophy flux equilibrium the result is S^iR) oc R 2 . The 
basis for further development is the immediate observa- 
tion that for d = | the scaling exponent for energy flux 
equilibrium, x e = 2, coincides with the the exponent 2 
of the equi-partition of enstrophy. Accordingly for d = | 
the law (j4|) is in agreement with enstrophy equipartition, 
and therefore also with a Gaussian pdf for the velocity 
differences. In the rest of this Letter we show that in 
dimensions | < d < 2, the flux remains small (for given 
total energy of the system) and the pdf 's do not change 
much from Gaussianity. 

A measure for deviations from Gaussianity is the skew- 
ness (|5|) which is now d-dependent, S(d). To compute it 
we need to separately find S 2 (R) and Sa(R). We start 
with the former. The structure function S 2 (R) can be 
computed from Eqs. (10) and (p^). In 2-dimensions 



in terms of the third order correlation Fi 



kqp • 



m = J 



C e e 2 ' 3 
2tt 2 



d 2 k 



exp(z/ci?cos< 



- I 



2 sin 2 4> k n e k 



sm 



dk 
p75 



[1 — exp (ikR cos cj>k)} 



= C e {eRY /i A 2 



o 

2/3, 



(21) 



An 



: <t>kd(f> 



2tt 2 



dn 



[1 — COs(KCOS</>fc)] 



27T(2/3) 
2 2 / 3 32 7rr(l/3) 

In d dimensions we write 



0.0855 



S 2 (R) = C 2 (d)(eR) 2/3 = C £ (d)(eR) 2/3 A 2 (d) . (22) 

However, it is easy to see that A 2 (d) does not become 
critical at d = |. On the other hand we will show 
that C e {d) (and therefore S 2 (R)) diverges when d — > |. 
Therefore we will estimate A 2 (d) by its value A 2 at 2- 
dimension, C 2 (d) = C e (d)A 2 . 

On the one hand S 3 {R) = 12eR/d(d+2) exactly. On 
the other hand the function Ss(R) is related to Fk qp by: 

a /d^ f d d kd d qd d p 

bz(R) = J — j^A2d~ o(fe + g + p)sin0feSin0 g sin0p 

xIm{[exp(ifci?cos</ifc) — l][exp(iqRcos(f) q ) — 1] 

x [exp(ipR cos cf)p) - l]}F kqp . (23) 

This is as far as we can proceed exactly. Now we will ex- 
press the third order correlator F kqp in terms of the sec- 
ond order n k . It is well known that this cannot be done 
without closure approximations. The latter are known to 
provide semi-quantitative estimates of the coefficients of 
correlation functions, and in the present context we ex- 
pect such approximations to perform better than in 3D 
due to the existence of the small parameter S(d) that we 
will expose in this calculation. The starting point is the 
equation of motion of n k (t) , which can be exactly written 



dn k /2dt 
h 



d d qd d p 
(27r) d 



5{k + q + p)V kqp F kqp . (24) 



A standard closure approximation expresses F kqp via n k . 
Proper closures in turbulence involve the following steps: 
first, one considers the sweeping- free renormalized per- 
turbation theory to first order. Second, one assumes a 
simple analytic form for the time-decay of correlation and 
response functions. A typical result reads |h], pj| 

F kqp = N kqp 9 kqp , (25) 
N kqp = V kqp n q n p + V pkq n k n q + V qpk n p n k . 

Here 9 kqp is a closure dependent "triad-decorrelation 
time" ; assuming simple exponential decay for the 2nd or- 
der correlation function and response function, one finds 



h qP = l/(7fc +lq + 7 P ) 



(26) 



In Eq. (26) j k is the width of the assumed Lorentzian 
line shape. The latter can be estimated as follows: 



lk = A 1 {d)ky / S 2 {l/k) = C 7 (d)e 1/3 fc 2/3 . (27) 

Here A 1 (d) is a coefficient of the order of unity. Since 
S 2 (R) diverges at d — | we will evaluate below A 1 {d) 
via its 2D value A 7 : 

C 7 (d) =Ayy/(M3j . (28) 



We should stress that although Eq. (g5|) is derived 
using usual uncontrolled closure approximations, once it 
is substituted into Eq. (|I|) the latter conserves the en- 
ergy and enstrophy invariants defined in Eq.(|l7j) in all 
d dimension. This distinguishes our analysis from some 
previous theories like Q which conserved enstrophy in 
2-dimensions only. Note that our equations of motion 
exhibit the equilibria (|l8| ) and (|lj]), ( po|) as exact results. 
Thermodynamic equilibria follow directly from the Ja- 
coby identities (|i~2| ) and (13) which yield N kqp = F kqp = 
and hence also I k — 0. To show that also the flux equilib- 
ria are satisfied exactly one needs to use the Kraichnan- 
Zakharov transformation |l2] ] and the Jacoby identities. 
There F kqp does not vanish in general. 

We now note that in dimension d = | when n k cx 1/fc 2 , 
F kqp vanishes by itself. This follows from the fact that 
N kqp vanishes due to the Jacoby identity (|i"3|). It is 
worthwhile to stress that in our context this result is 
derived only to first order in renormalized perturbation 
theory. It is however a stronger result which can be es- 
tablished order by order to all orders. 

Substituting (25) with n k — nf of Eq. ([l9|), and Eq. 
(f26|), we can rewrite Ss(R) in the form 



S 3 (R) = [C £ (d)] 2 Md)(£R)/C 7 (d) 



(29) 



where A 3 (d) is: 
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f d d kd d qd d p S(k + q + p)S^ § sin <£ fe sin <?!><? sin <^ r • ,t j, w • , % , • /- - yi 

Md) = J {2ir)2d - 2/3 - , 2/3 - ( ^-- )2d+5/3 [ Mk cos fo) + sm( g cos <t> q ) + sm(p cos <f> p )] 

X2 [fc d+2 / 3 (p 2 - q 2 ) + p rf+2 / 3 (g 2 - fc 2 ) + q d+2 /\k 2 - p 2 )] . (30) 



Combining Eqs. @, @, © and @ we find 



l S(d)=A 3 (d)/A 2 A 7 



(31) 



Obviously, A${d) = Q at d = |. This implies that 52 (i?) 
(i.e. C e (d)) diverge at this point. It is easy to prove that 
the first derivative of A3 (d) with respect to d at d = | is 
finite. We can therefore approximate ^(d) as 



A 3 (d) « a(d- 4/3) 



(32) 



up to orders of (d — |) 2 . On the other hand by direct 
numerical integration we found 

A 3 = A 3 (d)\ d=2 « 3.556 x 10~ 4 . (33) 

Thus in 2-dimensions we estimate 

S = 5(d) I d=2 « 0.0486/A 7 . (34) 

The experimental observation is that iS exp ~ 0.03. Tak- 
ing liberty to use this finding we estimate Ay w 1.62. 
This is in agreement with our expectation that A 1 is of 
the order of unity. Now we can estimate S(d) in the whole 
range | < d < 2 by using the linear approximation (32), 
and find 



0.0729 ^ 4 
<S(d) « — ^ — (d- 3) 



0.045(d- -) 
v 3 y 



(35) 



The main conclusion of this Letter is that although d = 2 
is finitely removed from d = | , the relevant small param- 
eter remains small all the way to d = 2, because of the 



numerical smallness of the ratio A3/A 2 , ps 0.0486. This 
smallness stems from generic geometric cancellations in 
the last line of the integrand in (|30|) . This originates from 
the structure of the vertex Vk qp and is therefore funda- 
mental to Euler equations in 2-dimensions. 

In summary, we have addressed the experimental find- 
ings of the statistics of 2D turbulence. It is hardly sur- 
prising that 2D turbulence is not strongly intermittent; 
in the inverse cascade regime there are no mechanisms 
to create rare events that are intimately related to the 
sharpening and enhancement of fluctuations as they are 
transferred down the scales in generic 3D direct cascades 
|i3|| . On the other hand the fact that the statistic of 2D 
turbulence are so close to Gaussian came as a major sur- 
prise. In this Letter we offered an explanation to this 
finding. We have identified d = | as a convenient point 
around which to develop a theory of 2D turbulence. The 
statistics there are Gaussian, but the spectrum in the in- 
verse energy flux regime is K41. The skewness is zero, al- 
lowing a sensible closure theory for d slightly larger than 
I . Estimating the skewness as a function of d, we are led 
to conclude that it remains small also in 2-dimcnsions. 
We can thus interpret 2D turbulence as a state close to 
equilibrium. In future work we will examine further the 
structure of the theory in the range | < d < 2 to assess 
further the quality of the closure approximation. 
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